OUTP-11-48P 

The fine-tuning of the generalised NMSSM 



Graham G. Ross, Kai Schmidt-Hoberg 

Rudolf Peierls Centre for Theoretical Physics, University of Oxford, 
1 Keble Road, Oxford 0X1 3NP, UK 

(N 

o 

(N 

O , Abstract 

We determine the degree of fine-tuning needed in a generalised version of the 
NMSSM that follows from an underlying Z4 or i?-symmetry. We find that it is 
significantly less than is found in the MSSM or NMSSM and extends the range of 
Higgs mass that have acceptable fine-tuning. Remarkably the minimal fine-tuning 
O { is achieved for Higgs masses of around 130 GeV. 
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1 Introduction 



Supersymmetric extensions of the standard model (SM) offer a solution to the hierarchy 
problem and can therefore generate a consistent Grand Unified or superstring unified 
model of the fundamental interactions. The remarkable success of the resultant predic- 
tion for gauge coupling unification has prompted an enormous body of work on SUSY 
phenomenology. The most immediate prediction is that there should be a spectrum 
of supersymmetric states with masses of the order of the electroweak breaking scale. 
However, none of the superpartners have been discovered yet and this has given rise to 
the little hierarchy problem, namely the difficulty in explaining why SUSY states should 
be so much heavier than the W and Z bosons. 

One may obtain a quantitative measure of the little hierarchy problem via a fine- 
tuning measure that quantifies the degree of cancellation between uncorrelated terms in 
the Lagrangian needed to achieve the mass splitting [l][2] . The degree of fine-tuning has 
been widely studied in the context of the minimal supersymmetric extension of the SM 
(MSSM). In the context of gravity mediation most of the low-fine-tuned areas of the 
model have already been ruled out by the non-observation of the SUSY partners and 
the recent LHC results leave only two significant regions to be checked. One of these is 
accessible to the LHC in its first 7 TeV run with lfb -1 of data while the second requires 
the full 14 TeV CM energy or improvement by one or two orders of magnitude sensitivity 
in the dark matter searches pi. For the case of gauge mediation the fine-tuning is worse 
and over all parameter space the fine-tuning required is greater than 1 part in 100 |4j. 

Given this, it is of particular interest to ask whether non-minimal SUSY exten- 
sions can ameliorate the fine-tuning needed and, if so, to determine the changes in 
phenomenology that may be implied. It has been known for some time [5] that the 
next-to-minimal SUSY extension that involves an additional SM singlet superfield, S, 
does reduce the fine-tuning, particularly at low tan/3, see also |6] for a recent analysis. In 
practice the studies have almost all been done assuming the theory has a non-i? symme- 
try that restricts the couplings in the superpotential to a trilinear singlet term together 
with a singlet coupling to the Higgs superfields (NMSSM); for general reviews of the 
NMSSM see [7|[8]. As a result the \i term is only generated when the scalar component 
of S acquires a vev, solving the \x problem as this term is of the visible sector SUSY 
breaking scale. 

Recently an elegant alternative solution of the /i problem has been proposed [9,10 



that uses a discrete i?-symmetry to forbid the bare fi term. A /i term of the correct 
magnitude is then generated by gravity mediation when SUSY and the i?-symmetry are 
broken in a hidden sector. Moreover the same symmetry forbids baryon- and lepton- 
number violating terms (apart from the Majorana mass term for right-handed neutrinos) 



1 



up to and including the non-renormalisable terms of dimension 5 in the Lagrangian - a 
significant improvement on the usual i?-parity of the MSSM. 

For the case that the MSSM is extended to include an additional singlet the same 
i?-symmetry allows additional terms in the superpotential compared to the NMSSM. 
Interestingly it has already been noted that the presence of such additional terms can 



further reduce the fine-tuning needed 11 13 . In this paper we extend these analy- 



ses to estimate the fine-tuning needed over the full parameter space available with the 
additional superpotential terms and find a significant reduction in the fine-tuning, par- 
ticularly for Higgs masses around 130 GeV. 



2 Singlet extensions of the MSSM 

The most general extension of the MSSM by a gauge singlet chiral superfield consistent 
with the SM gauge symmetry, the GNMSSM, has a superpotential of the form 

W = Wvukawa + + \S)H u H d + £S+ l -^ s S 2 + ^KS 3 (1) 

where WVukawa is the MSSM superpotential generating the SM Yukawa couplings. Al- 
lowing for SUSY breaking the soft terms associated with the Higgs and singlet sector 
are 

Kioft = m2 s \ s \ 2 + m lJ h u\ 2 + m 2 hd \h d \ 2 

+ h u h d + \A x sh u h d + ^kA k s 3 + ^b s s 2 + £ s s + h.c. J . (2) 

Here large letters refer to superfields while small letters refer to the corresponding scalar 
component. 

The free parameters of the GNMSSM relevant to the scalar sector at the electroweak 
scale are given by A, A\, k, A k , //, &//, fi s , b s , £, £ s and m\ u1 m 2 hd1 m 2 . We can trade the 
three soft masses for the Higgs vevs v, v s and tan f3. As discussed in Section [3j radiative 
corrections introduce a further dependence on the stop mass, M- t . Fixing v to give the 
observed gauge boson masses leave 13 independent parameters. 

2.1 The NMSSM 

The problem with the most general set of couplings appearing in eq. Q is that the 
gauge symmetry and supersymmetry do not constrain the dimensionful terms to be 
small. However to be viable they should be of the order of the supersymmetry breaking 
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scale or less; in the case of the /i term this is the "// problem" that also applies to the 
MSSM. A simple solution is to require that the renormalisable couplings are constrained 
by a I13 symmetry acting on the singlet chiral superfield and the Higgs superfields. Its 
superpotential has the restricted form 

Wnmssm = VVyukawa + \SH u H d + • (3) 



It is parameterised by a subset of the parameters of the GNMSSM given by A, A\, k, A 
v s , tan (5 and M^. In particular the bare mass terms are absent and an effective fi term 
is generated once the s field acquires a vev. Allowing for the soft SUSY breaking terms 
the latter is of order the SUSY breaking scale as is required for a viable theory and so 
the structure provides a nice solution to the fi problem. 

To date most studies have concentrated on this version of a singlet extension of 
the MSSM and we shall refer to it as the NMSSM. Of course, like the MSSM, an 
additional %2 -R-parity must be imposed to forbid the baryon- and lepton- number 
violating operators of dimension 3 and 4. Even with this there remain troublesome 
dimension 5 operators that, to inhibit proton decay, must be suppressed by a very large 
mass scale some 10 7 times the Planck mass! Also, to avoid problems with stabilisation 



of the hierarchy 14 and domain walls 15 an additional /^-symmetry must be imposed 
14lfT6l. 



2.2 The GNMSSM 



As mentioned above, an alternative solution to the /i problem has recently been de- 
veloped (9, 10 based on a discrete i?-symmetry. In the context of the general singlet 
extension of the MSSM the symmetry also explains why the other dimensionful terms 
are also constrained to be of order the SUSY breaking scale as is necessary for a viable 
model. A further advantage of these symmetries is that they also forbid all baryon- 
and lepton- number violating terms up to and including dimension 5 terms in the La- 
grangian. As a result there is no need to introduce a super-Planckian mass scale. Due 
to the i?-symmetry, there is also no need to enlarge the symmetry to avoid problems 



with the stability of the hierarchy and domain walls 10 



Two discrete i?-symmetries were identified, 1^ and Z^, under which the fields trans- 
form as given in Table IT Before SUSY/i?-symmetry breaking the superpotential is of 
the NMSSM form in eq.(3|. However after supersymmetry breaking in a hidden sector 
with gravity mediation dimensionful terms are generated. With these the renormalisable 
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Table 1: Charge assignments under Zjj^ for M = 4, 8. The labels 10 and 5 refer to the SU(5) 
matter content. 



terms in the superpotential take the form 10 



>V 2 fl rv W N MSSM + ™3/2 S + m 3/2 S 3 + m 3/2 # u #d , (4) 

W Z H ~ Wnmssm + "^3/2 S (5) 

where the ~ denotes that the dimensional terms are specified up to 0(1) coefficients. 

One may see that W Z fl is of the form of eq.Q, the GNMSSM superpotential, but 
with the constraint that the scale of the dimensionful couplings is set by WI3/2, the scale 
of supersymmetry breaking in the visible sector. After eliminating the linear term in S 
by a shift in its vev we obtain the general form of the superpotential given by 

Wqnmssm = >V Yu kawa + (a* + ^S)H u H d + ^ S S 2 + ^kS 3 (6) 

with /i ~ fi s ~ 0(m 3 / 2 ). We will refer to this model as the GNMSSM. In the case the 
underlying symmetry is 1^ we obtain the same form for the superpotential but with 
the constraint \xj\i s = \/(2k). 

In [l3] a slightly restricted form of the GNMSSM was studied. Their model, referred 
to as the SMSSM, has a restricted set of the parameters appearing in eqs. Q and @, 
corresponding to setting k and A K to zero. 



3 The scalar sector 

Fine-tuning refers to the difficulty in explaining why the electroweak breaking scale is 
so much less than the SUSY breaking scale. In order to provide a quantitative measure 
of fine-tuning it is necessary first to study the scalar potential that determines the elec- 
troweak breaking scale. In addition to the terms which originate from the superpotential 
and the soft terms given above (cf. eqs. (II]) and (|2j), we also have to take into account 
the .D-term potential which is given by 

Vn = (\K\> - \h£) 2 + 9 j\hlh u f ■ (7) 



4 



In our convention the Higgs fields are expanded around the neutral real vevs as 

h° u = v u + ^={Re(h° u ) + ilm(h° u )) (8) 
etc. such that 

v 2 = vl + v 2 d ~ (174 GeV) 2 . (9) 

It is well known that radiative corrections significantly affect the Higgs potential. For 
simplicity we only consider the leading one-loop contribution. Following [12] we add a 
term 



Uh u f = ^q^W (io) 



to the Higgs potential, with 

8 = In — ^ = f . In — ^ . 11 

Here M f ~ is the average stop mass while m t and h t are the top mass and Yukawa coupling 
respectively. 

The condition for the scalar potential to develop a minimum at v u ,Vd,v s can be 
written as 

m 2 hu + (Xv s + fi) 2 - \v 2 {g 2 - 4A 2 ) cos 2 (3 

- (bfi + X(v s (A x + v s k + /i s ) + 0) cot f3 + ^g 2 v 2 (1 + 6) sin 2 (3 = 0, (12) 
m\ d + (A Vs + /i) 2 - ^ V - 4A 2 ) sin 2 (3 

- (bfj, + \{v s (A x + v s k + fi s ) + £)) tan /3 + -#V cos 2 /5 = , (13) 

(f 2 A/i + /i s £ — £ s v 2 \(A\ + // g ) cos (3 sin /3) + (b s + m 2 + t> 2 A 2 + ^ 2 

+ 2k£ - w 2 kA sin(2/3))t; s + + 3/i s )w 2 + 2k\ 3 = . (14) 
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The mass matrix of the CP even Higgs fields in the basis (h u , h^, s) can be written 

as 

Mil = cot{/3){\{v s {A x + v s k + fjL a ) + £) + tyx) + Mfsin 2 (/3)(1 + 5) , 
M 2 2 = tim(f3)(\(v s (A x + v s k + fi 8 ) + + bfi) + Mj cos 2 (/3) , 

M| 3 = v,k(A k + 4u s « + 3/x s ) + — (v 2 X(A x + p,) sin(2/3) - 2v 2 \fi - 2^ - 2£ s ) , 

M 2 2 = -\{v s {Ax + v,k + H.) +£)-b(i+ {v 2 X 2 - M 2 /2) sin(2/3) , 
M 2 3 = fA(2sin(/3)(w s A + //) - cos(/3)(A A + 2v s k + /x s )) , 
M 2 3 = t>A(2(> s A + /i) cos(/3) - (A x + 2kw s + /x a ) sin(/3)) , 

where the soft masses m\ u ,m 2 hd and m 2 have been eliminated through the minimisation 
conditions. Sometimes it is useful to consider the basis in which only one of the (doublet) 
Higgs h obtains a vev. This is done by a rotation of angle (3 between the two doublet 
Higgs fields. In this basis (h, H, s) the mass matrix then reads 

M\ x = Mf(cos 2 (2/3) + 5sin 4 /3) + AV sin 2 (2/3) , 

K _ = b±± m_M^^ + + j, _ AV) M , 

M 3 2 3 = v s k(A k + Av s k + 3// s ) + — (VA(A A + /i s ) sin(2/3) - 2w 2 A / u - 2/i s £ - 2£ s ) , 

M 2 2 = -|(M| - AV) sin(4/3) + 5M 2 Z cos /5 sin 3 /3 , 
M 2 3 = wA(2w s A + 2/i - sin(2/3)(A A + 2w s k + /i s )) , 
M 2 3 = -wAcos(2/3)(A A + 2v s k + fi s ) . 

The smallest eigenvalue is always bounded from above by the smallest diagonal entry. In 
particular M 2 X corresponds to the well-known upper bound on the lightest Higgs mass 
in the singlet extended MSSM. In comparison to the MSSM case there are additional 
terms which depend on the trilinear coupling A. One can see from the form of M\ x that 
the lightest Higgs mass can become larger for larger coupling A. However, the coupling 
A is bounded from above by the requirement of perturbativity up to the GUT scale. 
For tan/3 > 2 this roughly corresponds to A < 0.7, cf. e.g. pj. Furthermore, a large 
value of A generically leads to sizable mixing between the singlet and doublet Higgs 
fields, which decreases the lightest Higgs mass. This can be avoided by tuning Ax such 
that the matrix element M\s and hence this mixing becomes zero. Note, however, that 
starting from a specific high energy theory that determines the boundary values of the 
parameters it may not be possible to do this. 
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Once we fix tan /3 the elements M\ x and M\ 2 cannot be varied independently. How- 
ever the remaining five independent elements can be varied independently by adjusting 
the free parameters. The situation is the same in the usual NMSSM and for its gen- 
eralised version. In this sense the new parameters of the generalised NMSSM don't 
introduce additional freedom in the mass matrix. Nevertheless the corresponding fine- 
tuning can be different in both cases and will be studied in the following. 



4 Fine- Tuning 



As introduced in [T}[2], a quantitative estimate of the the fine-tuning of the EW scale 
with respect to a set of independent parameters, p, is given by 

A = maxAbs[AJ, A p = = • (15) 

omp v 2 op 

The quantity A -1 gives a measure of the accuracy to which independent parameters must 
be tuned to get the correct electroweak breaking scale. The independent parameters are 
usually taken to be those at a high energy scale determined by some underlying theory 
relevant at that scale. To determine them at the electroweak scale one has to compute 
their renormalisation group running. 

While the form of the GNMSSM superpotential is expected quite generically, the 
soft supersymmetry breaking terms as well as the associated fine-tuning depend on the 
details of the supersymmetry breaking sector. In what follows we choose a simple scheme 
for supersymmetry breaking, similar to the well-known CMSSM. We take the high scale 
parameters to be a common gaugino mass, raj/2, a common A-term, A , bilinear and soft 
bilinear terms for the doublet and singlet Higgs fields, /io, &/io, p>®, b®, and a common 
scalar mass mo for all scalar fields. The parameter set is completed by the Yukawa 
couplings Ao and KoQ We therefore have 9 free parameters in the GNMSSM case. To 
obtain the relations between electroweak and the high scale parameters we solve the 
one-loop RGEs |8] numerically for different GUT scale boundary conditions and use it 
to interpolate over the full parameter space. The corresponding numerical relations we 
use can be found in the Appendix. In the case of the NMSSM it is well known that 
with universal mo, m\i% and Aq the Higgs mass cannot be enhanced with respect to the 
MSSM case and we will therefore not consider the NMSSM case any further. A more 
general structure of the SUSY breaking terms which also allows an enhancement of the 
Higgs mass in the NMSSM will be presented elsewhere. 



1 Using the modified definition for fine-tuning [17], appropriate for measured parameters, the contri- 



bution from the top Yukawa coupling is sub-dominant 18 
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4.1 Fine-tuning scan 

At the low scale three input parameters may be traded for v, v s and tan (5 using the 
minimisation conditions. We choose these parameters to be /xo, b/io, b® for the GNMSSM 
and Ho, bfio for the MSSM. We performed a scan over the remaining parameter space of 
the MSSM and the GNMSSM to compare their corresponding fine-tuning. The result 
of these scans is shown in Figure [TJ 




100 110 120 130 140 100 110 120 130 140 



m hl [GeV] m hl [GeV] 

Figure 1: Comparison of fine-tuning in the MSSM (orange) and GNMSSM (blue) for a minimal 
unified gaugino mass m\/2 > 300(500) GeV. For the GNMSSM the lowest fine-tuning is at 
low tan/3 = 0(2) while in the MSSM it is at higher tan/3 = 0(10). Remarkably the lowest 
fine-tuning is achieved for rather large Higgs masses, ~ 130 GeV. 

In both cases we scanned over 300 (500) < m 1/2 < 1000 GeV, < m < 1000 GeV, 
< A < 1000 GeV, and 2 < tan/3 < 10. In addition, in the GNMSSM we scanned 
over 0.5 < A < 2, < k < 0.5, -1000 GeV < v s < 1000 GeV, and -5000 GeV < /jP s < 
5000 GeV. Note that we restricted the effective /i-term to /i e ff = (/i + Xv s ) > 100 GeV 
to ensure that the charginos are above the experimental limit. The other superparticle 
masses are predominantly determined by the universal gaugino mass mi/2- As a rough 
estimate for the current bounds from the LHC we impose a lower limit of > 300 GeV 
and to demonstrate how things will change if the bounds improve we also show the case 
m i/2 > 500 GeV. A more detailed numerical study of these constraints will be presented 
elsewhere. 

It may be seen that the fine-tuning in the GNMSSM is reduced relative to the MSSM 
particularly at larger Higgs masses. Remarkably, the lowest fine-tuning is achieved 
for Higgs masses of ~ 130 GeV. Comparing the cases for = 300 GeV and 
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m i/2 — 500 GeV it may be seen that although the overall fine-tuning increases with 
mi/2, the Higgs mass corresponding to the smallest fine-tuning is almost unchanged. 

With the given universal soft parameters at the GUT scale we don't find any viable 
points with large Higgs masses and small values of the singlet mass parameter fi®. This 
can also be seen in Figure [2j where we show the lightest Higgs mass versus pP s for points 
which have a fine-tuning below 100. This means that the corresponding singlet states 




90 100 110 120 130 



m h] [GeV] 

Figure 2: The lightest Higgs mass versus the high-scale singlet mass parameter for points 
with fine-tuning A < 100. 

will be rather heavy leading to a phenomenology very similar to the MSSM, but with 
a smaller fine-tuning and a heavier Higgs than would be expected from the masses of 
the superpartners. It also implies that the slightly more 'NMSSM-like' version following 
from an underlying 1^ symmetry will have somewhat larger fine-tuning than the more 
general 1^ case, because the relation \xj ' \x s = A/ (2 k) implies a fi term which is of similar 
size as /z s . Note that both these conclusions crucially depend on the assumptions we 
have made about the supersymmetry breaking terms and could be very different in a 
more general setup. In fact it is known that in the NMSSM with more general SUSY 
breaking terms the singlet states can be light leading to interesting implications for 
phenomenology. We will consider these more general SUSY breaking terms for the 
GNMSSM in a future work. 

4.1.1 Origin of the reduction in fine-tuning 

There are two main reasons for this reduction in fine-tuning. The first, which applies to 
both the NMSSM and GNMSSM, is that there is a contribution, AV sin 2 (2/3), to the 
matrix element, M\ x , that sets the upper bound on the lightest Higgs mass. This has 
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the general effect of shifting the fine-tuning envelope in Figure [T] to the right; how much 
depends on whether it is possible to saturate the bound with parameters maintaining 
low-fine-tuning. Due to its additional parameters this is easier in the GNMSSM. Note 
that the effect is significant only for low tan (3 (tan /3 ~ 2) . 

The second, which applies mainly to the GNMSSM, is the additional stabilising 
term in the Higgs potential that the additional singlet coupling generates and reduces 
the relative magnitude of the Higgs vevs to the SUSY breaking scale. This is most 
clearly seen in the large /x s limit where one can integrate out the S superfield at the 
supersymmetric level; this limit is applicable only to the GNMSSM and not to the 
NMSSM. In this limit one obtains the superpotential term \ 2 (H u H<i) 2 / '/x s . The leading 
contribution to the Higgs potential coming from this term is given by 

V ~ 4\ 2 (\v s + fi)h u h d (\h u \ 2 + \h d \ 2 )/fx s . (16) 

This is proportional to sin 2j3 and so is peaked at low tan (3 although not as sharply 
as for the additional Higgs mass contribution. The presence of this term was found [12] 
to significantly reduce the fine-tuning relative to the NMSSM and explains the origin of 
the difference seen in Figure [T] for the case /x s is large. For lower values of /j, s the singlet 
coupling generates stabilising terms in both the NMSSM and the GNMSSM although 
they are likely to be more important in the GNMSSM case. 



5 Summary and Conclusions 

Given the enthusiasm for SUSY as a solution to the hierarchy problem there has been 
surprisingly little attention paid to the general version of the NMSSM. This has been 
due to the apparent un-naturalness of the model which has additional mass terms in 
the Lagrangian that must be no larger than the SUSY breaking scale in the visible 
sector. However the recent discovery that simple discrete i?-symmetries can naturally 
achieve this has made the GNMSSM as natural as the NMSSM. Indeed the fact that 
the symmetry also eliminates the dangerous dimension three, four and five baryon- and 
lepton-number violating terms in the Lagrangian and avoids destabilising tadpoles and 
domain wall problems renders it a more promising starting point than the NMSSM. 

In this paper we have determined the fine-tuning in the GNMSSM for the case 
of gravity mediated SUSY breaking. We found that the fine-tuning is significantly 
reduced compared to the MSSM and in particular significantly extends the range of 
Higgs mass that have acceptable fine-tuning. Interestingly the region with smallest 
fine-tuning corresponds to Higgs masses of ~ 130 GeV. 
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For the universal boundary conditions we consider the singlet mass scale entering in 
the superpotential was typically rather large compared to that of the other dimensionful 
parameters, implying that the phenomenology is very similar to that of the MSSM, 
although the fine-tuning can still be significantly reduced. 

Our analysis of the GNMSSM has been somewhat rudimentary, relying on analytic 
approximations for the RGE running of parameters. Work on a more complete numerical 
analysis is in the planning stage and we hope to be able to present a more detailed study 
of the fine-tuning and associated collider phenomenology in the future. 
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A RGE running 

The independent high-scale parameters are given by m 1 / 2 , A , /x , bfj,o, l^,b®, mo, m^m^, 

A and kq. Using the one-loop RGEs the relation between electroweak and GUT scale 
parameters can, for the low tan/3 region of interest, be approximately written as 

m 2 hu ~ -0. 10 Al - 0. 77 m 2 + 0.39 A m 1/2 - 2.57 m 2 1/2 + 0.58 (m° h J 2 

+A [0.21 m{ /2 - 0.08 (m° h f] (17) 

m\ d ~ 0.44m 2 /2 + 0.94 (m°J 2 - A [0.07m 2 + 0.11 (m°J 2 + 0.10 (m°J 2 ] (18) 

H ~ 0.97/io -0.18Ao/x (19) 

bfi ~ 0.93 bfi - 0.40 A fi + 0.22 m 1/2 fi + \ [0.07 A pi - 0.38 b/i } (20) 
m 2 s ~ 0.90m 2 + A [0.23 KJ 2 - 0.19 KJ 2 - 0.10m 2 ] 

-0.45/t m 2 + 0.14A /t m 2 (21) 

H, ~ 0.97^° - 0.36 A /i°- 0.30 k fi° 8 + 0.12 A kq fj? s (22) 
b s ~ 0.93 b° s - 0.11 A fi° s - \ [0.298 6° + 0. 10 A /i°] 

-k [0.38 b° s + 0.05 A n° s ] + A k [0. 17 b° s - 0.07 bfi ] (23) 

A ~ 0.78 A -0.21 A 2 (24) 

k ~ 0.58 K -0.16^ (25) 
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m 



A\ ~ 0.54 A + 0.23 m 1/2 - 0.31 \ A - 0.16 k A (26) 

A K ~ 0.84A + A [0.10 m V 2- 0.52 A ] -0.42 k A + 0.21 A k A (27) 

~ 0.74m 2 + 0.14A mi /2 + 5.24m 2 /2 -0.14(m£j 2 + 0.08A m 2 /2 (28) 

m 2 tR ~ 0.53 m 2 + 0.25 A m 1/2 + 3.90 m 2 /2 - 0.20 (m£j 2 + 0. 13 A m 2 /2 (29) 

A t ~ 0.24 A - 2.14 m 1/2 - 0.12 A m 1/2 (30) 
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